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ABSTRACT 

A visual problem-solving technique a^>plicable to 
several different classes of mechanics time-dependent problems is 
discussed. The computer is used to solve the equations of motion of 
various mechanical systems by one of several standard methods^ and 
the solutions are displayed in time increments. A specific example is 
provided to illustrate this technique. (MLH) 
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Pictorial Solutions in Advanced Mechanics 

Abstract 



Q)ue to the , expansive amount of lecture theory workload 
and the lack of time available for peripheral 'probl em^-sqlving/ 
innovative teaching techniques in advanced level mechanics 
courses frequently go unconsidered. Another difficulty is 
that at this level, few texts are standard' and even fewer 
instructors fo'llow these with exactness. The one aspect of 
advanced courses , in particular in advanced dynamics, which 
remains common to most classrooms is that the same problems 
are solved, or at least the governing equations are solved^ 
Due. to the complexity of such equations, however, a simplir- 
fied case" or two may be investig.ated briefly before new . 
theory is begun . • ' ' " 

What is described in this paper is a visual problem- 
solving technique applicable to several different classes of 
mechanics tirae-dependerft problems. The computer is used to 
solve the equations of motion of various mechanical systems 
by one of several standard methods and the solutions are 
displayed in time increments. With this technique"., impor- 
tant mechanical/physical parame-ters can be varied to give ^ 
parametric investigation .capabi lity to the classroom. 

The problem demonstrated here is a two mass coupled 
system whose aquations of motion form a s"et of 
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nonlinear, coupled dTf f erenti al equations of two degrees of 
freedom.- The solution to this problem has been studied^in 
detail, but mainly for cases where the nonlinear behavior is^ 
neglected, 0t4rerwise it generally defies classroom attention. 
First the initial conditions aVe prescribed and then the' 
Runge-Kutta method is applied to extend the solution for each 
time step. The generalized coordinate are stored' in s( com- 
puter and fin&lTy a graphtcs routine is invoked so that at 
uniformly s||)aced time steps a picture is displayed or printed. 
If an interactive computer facility is available, the results 
may be viewed directly then electronically redrawn as time is 
augmented, otherwise motion picture film or videotape mu.st be 
^ used to first copy each time frame and then to display the 

, complete motion. • * " * * - 

This technique can be useful in exploring the governing 
differential equations of many problems, in the area of 
dynamics and vibrations. In particular the motipn of other 
systems such as pendulums, gyroscopes, and rigid bodies is 
easily introduced into the classroom. All routines are 
written' in f^ortran and available upon . request . 
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Purpose 

The use of laboratory experiments to furnish a physical 
foundation for theory at the beginning level in Meckanics is 
well" establ ished. In fact it is customary to devote-formal 
course time nearly equally between lecture and laboratory, 
PractlcaTh^de-emphasis of laboratory however is eventually 
assumed by almost all students because cowrse credit is 
applied to that aspect of the educational process very spar- 
ingly. The purpose of this paper therefore is to suggest a 
general method whereby continuing physical emphasis can be 
maintained in the advanced level mechanics curriculum without 
an inordinate penalty in instructional preparation time. This 
purpose is demonstrated by examples of several solutions fol- 
lowing a presentation of the appropriate theory. It is hoped 
that the ease by which pictorial solutions may be applied to 
mechanics will generate, an interest on the part of educators 
in using this technique and in encouraging the student to 
devote at least some attention to this method of investigating 
advanced mechanics-* 

Theory 

The equations of motion for mechanical systems can be ' 
written generally as a set of dif f eren-tial equations in N 

i 

generalized coordinates, q . , as the dependent variables for 
systems with N degrees of freedom, i - 1,2, N, with the 
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one independent variable time. These systems are always 2nd 
order ordinary differential equations with varia^ble coeffi- 
cients, and non-1 inear in nature. Using the Lagrangian form 
ulation for.deriving the equations of Motion; it is seen that 
the Lagrangian function, L^is at worst a function of q-|, q-f, 
and t. Since single differentiation af the. Lagrangian with 
respect to time gives the highest order derivtttives of the^ 
q-j and q^'s (q-,-)* then eventually a syste^ of equations 
results which is linear in the q\-^s, or 

since _d 
dt 

then at worst 

A^-j(qj., q^.t) q'i = B.(q., q^- , t); * 

J 1, •••5 N# 

(paired subscripts denote summation) 
The A^-j coefficient matrix is not generally singular so that 

Ajj exists leading to a set of N equations for the q-j«s of 
tjfe form: 

When this form is established, an effective method 
must be utilized which by^uadrature generates the q^-'s as a 
function of time, A very general accurate and speedy process 
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is advantageous. Also, practical experience dictates "^he 
further condition that the ititegrati^gn step should be capable 
of being varied easily. These conditions fit the Runge- 
Kutta^ Method attributes probably better than any other. 

The advantage* of this method technically over others for 
the programmer, a peripheral pursuit for most of us, is that 
prior solution values q-jCtp), qj(tp) are sufficient to deter- 
mine^ the values qj(tp+i), q^Ctp+i). Thus values, for 
instance, at tp^i , *n-2> *n-3> unnecessary and 

need not complicate the integration ^process or clutter machine 
memory, an important consideration in core-bound computers. 

Although this method is normally stated for first-order 
differential equations, it may be ext'ended simply [3] by 
restating the original problem slightly. For instance if 
there are N q-j's and N equations with (2N) initial conditions 
qi(to)> ^i(*o)» t\\Br\ through the governing differential equa- 
tions in q-j , the solution q-jCt-j) and-q^*(ti) are determined 
from the solution at to with t] = to+ h. In the general 
Runge-Kutta Method [a third order example is used], a change 
of variable is employed to generate 2N dependent variables as 

qi ' 

. l^i = fi(qi' Uj, t); j = 1, N . 



with initial conditions ..Ui'tto)^ qi(to) ' 
The solution to this system at t] is: 

■ . r ^ , ' 

let ki = hUi(to) ki- = hFi-(to, qj(to), Uj(to)); J = N 
ly = hlUj^to) + ^ki-], SLi = hFi(to + Hh, qj(to)+ %kj , 
U.(to) + %k.); j = 1 N 

Mi = h[Ui(to) + |?i],Mi = hFi(t^ |H.qj(to).+ f Xj, 

then ^^{t■^) = q^ttg) + |. (2k^. + 31^ + 4m^. ) 

Ui(ti) = Ui(to) + I (2ki + 3£i + 4mi) 
at this point h = t2 - t-| , can be changed, and the process 

» 

repeated. 

Procedure 

This method is^yery useful if a routine exists which can 
transform the q-j's into a graphical solution at selected 
time steps. These time steps are chosen so tfiat^ ihey match 
framing speeds of conventional cameras, where typically h^.05 
sTjBc, or CRT or video equipment. In, some complex cases, the 
h used for the video work may be a multiple of the numerical 
solution time step. 

8 
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The computer coding subroutines exist at this writing for 
picturing circles, or spheres:, and' linkages or bars. These 
are adequate to gener.ate most systems commonly used in class- 
room examples- The main computer program solves the equations 
of motion for the number of time steps until another print 
step is reached then^he q^- 's are used to evaluate the center* 
point coordinates and radius of any circles to be, printed as 
well as the two end pin connections and width*of each linkage. 
Next, the main computer program calls these generation sub- 
routines' which clear the previous picture. from storage and 
form the new image. At this writing printer output has been 
utilized on a 10 inch by 10 inch grid with a caption, but the 
method is more efficient on CRT dispflays or on plasma screens 
where available. 

The computer program is available as written in FORTRAN 

IV, and requires only approximately 100 K bytes of core. Sug- 

O 

gested types of mechanical systems which can be easily 



depicted are: 
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Pendulums 
Tops 

N-body Central Force Problems 
Spring Mass Systems, 
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Sample Problem 

A classical problem involving two- masses coupled by a , 
string is consider^. One mass, m, is allowed to move on a 
frictionless table iilLhil«e a second mass, km, is attached 

through a hole in the table by means of an inextensible cord 

\ 




Two Mass System 

The equationS^'oif motion [4] are given by 

(1 + k) r - r9^ = - kg and Zr^'r + r*- 9 = 0 ■ and sub- 
ject to 9(to), 9(to), r(to), r(to). 
Letting the Runge-Kutta variables be 

q-| = 9, q2 r, = 9, = f, and N = 2; 
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The algorithm becomes: 

ki = h.lli(to) 

k2 = h U2{to) 

I] - h[Ui(to) + % ki;] 

■^2 = h[U2(to) h k2] 

mi = h[Ui(to) + ^hl ■ 

m2 = hCU2Cto) + I Izl 

ki = - 2hUi(to) U2(to)/q2(to) 
2 

. . . k2 = h(q2(to) (to) - kg)/(l + k)' 

£l 2h(Ui(to) + h ki)(U2(to) + k k2)/(q2(to) + k .^2) 
^2 = . h(Cq2(to)+ % k2)(UiCto) + % k])' - kg))/(l + k) 
mi = -2h(CUiCto) +|^T)(U2Cto) + J ^2) )/ (^2 CtoKt | ^2.) 
mz = .h(Cq2(to) + | •^2)(U Cto) + 3^£i)=^ - kg)/jl + k) 
q^it■^) = q^(to) + I {Zk^ + 3l■^ + ) 

qi{ti) = q2(t6) + 1 (2k2 + +-4m2) 




Ui(tT) = Ui(to) + I |2fei + 3£l + 4m]) 



- ir - 

Initial conditions used for the computer generated 
pictures were: 



rn = 1*0 feet 



= 10.0 feet/sec. 



Oo = 0.0 radians 

9q = 3-0 radians/sec 



4-0 



g = 32.2 feet/sec- 



r 



to (time at stai^t^f motion) = O.D sec. 
*max (tTJ^s at end of motion) = 10.0 sej6. 
h - •OO'l seconds 

Print step at every .05 seconds. 
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